Can Nonequilibrium Spin Hall Accumulation be Induced in Ballistic Nanostructures? 
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We demonstrate that flow of longitudinal unpolarized current through a ballistic two-dimensional 
electron gas with Rashba spin-orbit coupling will induce nonequilibrium spin accumulation which 
has opposite sign for the two lateral edges and it is, therefore, the principal observable signature 
of the spin Hall effect in two-probe semiconductor nanostructures. The magnitude of its out-of- 
plane component is gradually diminished by static disorder, while it can be enhanced by an in-plane 
transverse magnetic field. Moreover, our prediction of the longitudinal component of the spin 
Hall accumulation, which is insensitive to the reversal of the bias voltage, offers a smoking gun to 
differentiate experimentally between the extrinsic, intrinsic, and mesoscopic spin Hall mechanisms. 

PACS numbers: 72.25.Dc, 73.23.-b, 85.75.Nn 
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Introduction — When electric current flows along a con- 
ductor subjected to a perpendicular magnetic field, the 
Lorenz force deflects the charge carriers creating a trans- 
verse Hall voltage between the lateral edges of the sam- 
ple. The normal Hall effect is one of the most familiar 
phenomena, as well as a widely utilized tool, in condensed 
matter physics In the absence of external magnetic 
fleld, more esoteric Hall-type effects involving electron 
spin become possible in paramagnetic systems with spin- 
orbit (SO) couplings — the opposite spins can be sepa- 
rated and then accumulated on the lateral edges when 
they are transported by a pure (i.e., not accompanied 
by any net charge current) spin Hall current flowing in 
the transverse direction in response to unpolarized charge 
current in the longitudinal direction. For instance, the 
SO dependent scattering off impurities, which deflects 
spin-f and spin-| electrons of an unpolarized beam in 
opposite directions, and it is partially responsible for the 
anomalous Hall effect in ferromagnetic metals JJ, has 
been invoked in early studies to predict the extrinsic (i.e., 
due to impurity scattering) spin Hall effect 0- 

The pursuit of all-electrical spin current induction and 
manipulation in semiconductor spintronics Q has rekin- 
dled interest in the realm of the spin Hall effect (s). In 
particular, recent theoretical arguments have unearthed 
possibility for pure transverse spin Hall current that is 
several order of magnitude greater than in the case of 
the extrinsic effect, arising due to intrinsic mechanisms 
related to the spin-split band structure in SO coupled 
bulk or mesoscopic 0,00 semiconductor systems. 
While these theories are formulated in terms of not di- 
rectly observable spin currents (which are the Fermi-sea 
quantity and not conserved in infinite homogeneous sys- 
tems of the intrinsic effect 0, 0) or the Fermi-surface 
property and conserved ones in the transverse electrodes 
attached to mesoscopic samples 0,00), their detection 
requires to measure nonequilibrium spin accumulation 
that they would deposit at the sample edges I2L Ifll llfl . 
The controversy in theoretical interpretations [ll], of 
recent breakthrough observations of the spin Hall accu- 
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FIG. 1: (color online), (a) The out-of-plane component 
(Szir)) of the nonequilibrium spin accumulation induced by 
nonlinear quantum transport of unpolarized charge current 
injected from the left lead into a two-terminal clean 2DEG 
(of size L — 30a > Lso, ffl — 3 nm) nanostructure with 
the Rashba SO coupling tso = O.lto and spin precession 
length Lso ~ 15.7a. Panel (b) shows how lateral spin-t 
and spin-| densities will flow in opposite directions through 
the attached transverse ideal (tso ~ 0) leads to generate a 
linear response spin Hall current [ly]^ out of four-terminal 
2DEG (L = 8a < Lso) nanostructures H], which changes 
sign [IyY{—V) = —[IyY{V) upon reversing the bias voltage. 



mulation |^ is largely due to the fact that no theory 
exists that demonstrates the existence of such accumula- 
tion in ballistic SO coupled finite-size structures accessi- 
ble to experiments |9|, [l^l ■ 

Thus, to answer the question posed in the title 
of this Letter, we formulate a Landauer-Keldysh ap- 
proach [13I to the spin Hall accumulation problem in 
two-terminal devices, which treats phase-coherent spin- 
charge transport while taking into account all bound- 
aries, interfaces, and electrodes of the semiconductor 
nanostructure [l5| . We predict a nonequilibrium spin ac- 
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cumulation induced at the lateral edges of a ballistic two- 
dimensional electron gas (2DEG) with Rashba SO cou- 
pling attached to two ideal (interaction-free) semi-infinite 
leads as a response to unpolarized charge current flowing 
through the sample, as shown in Fig. ^a). The electric 
field along the z-axis, which confines electrons within the 
xy-plane of quantum well in semiconductor heterostruc- 
ture, generates an effective momentum-dependent mag- 
netic field Br{ (p) [which does not break time-reversal in- 
variance] due to the Rashba SO coupling j3, 0| • The 
pattern of local spin density throughout the device is ob- 
tained from 



(S(r)) = - 



^Tr,pi„[6-G<(r,r;£;,F)], (1) 



EF~eVl2 

where the exact spin-dependent lesser nonequilibrium 
Green function G<(r, r'; i?, is evaluated in the 
steady-state quantum transport through the semiconduc- 
tor nanostructure attached to external probes 0| . Here 
V is the applied bias voltage between the leads and ^a- 
is the spin-i operator. In the nonlinear phase-coherent 
transport regime, Fig.^a) and Fig. |21 demonstrate that 
transverse profile of the out-of-plane (S'z(r)) component 
of the spin accumulation develops two peaks of oppo- 
site sign at the lateral edges of the 2DEG. Upon re- 
versing the bias voltage, the edge peaks flip their sign 
{S^ir))_y = -{S,{r))y. Around the left-lead/2DEG 
interface, (5*2 (r)) ^ in Fig. ^a) is suppressed since 
unpolarized electrons arc injected at this contact. 

All three features of {Sz{v)), delineated here for phase- 
coherent transport regime at low-temperatures, are re- 
markably similar to the general phenomenology of the 
spin Hall effect demonstrated convincingly in two very re- 
cent experiments lioj . These experiments focus on the 
optical detection of the spin Hall accumulation of oppo- 
site sign on the lateral edges of two-probe semiconductor 
structures (in the semiclassical transport regime at finite- 
temperatures) : (i) weakly SO coupled 3D films of GaAs 
or strained InGaAs |9|; and (u) strongly SO coupled 2D 
hole gases The former experiment has been inter- 

preted as the manifestation of the extrinsic effect 0, 
due to small resulting polarization (~ 0.03%) and no SO 
splitting of the band structure (9| , while the later is con- 
sistent with the intrinsic mechanisms |^ 0| because of 
spin-split quasiparticle energies and much larger spin po- 
larization (- 1%) ml. 

However, it has been estimated that spin Hall accumu- 
lation induced by the extrinsic effect is far below present 
experimental sensitivity [l^ . On the other hand, it has 
been argued [ill that the intrinsic spin Hall current 
in the bulk, determined only by the equilibrium distribu- 
tion function and spin-split Bloch band structure 0, 0| , 
does not really transport spins, so that no spin accumu- 
lation at the edges is possible in the absence of impuri- 
ties • We bring such debates to an end in Fig. ^b) 
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FIG. 2: (color online) . The one-dimensional transverse spatial 
profile of the spin accumulation {S{x = 78a,y)) across the 
200a X 200a 2DEG with the Rashba SO coupling tso = 0.02to 
through which ballistic quantum transport takes place in the 
nonlinear regime eV = 0.4to (sohd lines) or the linear regime 
eV — W~^to (dotted hues). The width of the edge peaks of 
{Sz{x = 78a, y)) is ^ Lso/2 = Tvato/itso- 



by attaching two additional transverse ideal semi-infinite 
leads at the lateral edges of a perfectly clean 2DEG and 
showing how spin-t and spin-| densities will flow through 
those leads in opposite transverse directions to g enerate 
the spin Hall current driven by mechanisms 17] on the 
mesoscopic scale .18i] set by the spin precession length 
Thus, Fig.^b) demonstrates convincingly that spin Hall 
effect in SO coupled ballistic nanostructures can be used 
as an all-electrical semiconductor-based spin injector 0. 

The attempts |23| to understand the spin Hall ac- 
cumulation in macroscopic disordered Rashba spin-split 
2DEG attached to two massive electrodes have reached 
contradictory conclusions (vanishing (^^(r)) on the lat- 
eral edges lldl v s. {Sz{r)) 7^ at the edges and within 
the sample |20j), which can be traced to an incom- 
plete treatment of experimentally relevant measuring ge- 
ometry within the semiclassical diffusion equation ap- 
proaches. A plethora of spin transport phenomena il- 
lustrate the need to treat the whole device geometry 
due to the presence of SO couplings, even in the semi- 
classical transport regime |l8l| . For example, spin re- 
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laxation in confined Rashba structures is quite differ- 
ent from the bulk D'yakonov-Perel' (DP) spin relax- 
ation [3 because of the transverse confinement effects 
or chaotic vs. integrable boundary scattering l^^. Also, 
Br,(p) in Rashba SO coupled wires is almost parallel 
to the transverse direction (in contrast to the infi- 
nite 2DEG where no unique spin quantization axis ex- 
ists 0). Mesoscopic transport techniques, developed to 
treat the whole measurin g ge ometry as demanded by 
quantum coherence effects are well-suited to handle 
all relevant details of the spin Hall transport measure- 
ment setups 1^ Hot . Therefore, we apply the Keldysh 
Green functions 13] to the Landauer-type geometry 
where finite-size sample is attached to two macroscopic 
reservoirs via semi-infinite ideal leads (which simplify 
the boundary conditions) — here the current is limited by 
quantum transmission through a potential profile while 
power is dissipated non-locally in the reservoirs 0, 0| . 

Landauer- Keldysh Green function approach to 
nonequilibrium spin accumulation. — The effective mass 
Hamiltonian modeling the ballistic finite-size 2DEG 
with Rashba SO coupling in Fig. ^ is given by 0,0] 

p^ -\- (X 

H = '^^^^ ^ + {pyCTx - Px^y) + l^conf (a;, y), (2) 

where (ct^;, fJj,, cTz) is the vector of the Pauli matrices, 
(px,Py) is the momentum operator in 2D space, a is the 
strength of the Rashba SO coupling 07 and Vconf{x,y) 
is the transverse confining potential. In order to evaluate 
the nonequilibrium Green functions for a sample of arbi- 
trary shape attached to the ideal leads, we follow Ref. [l^ 
and employ the local orbital basis. In this representation, 
the Rashba Hamiltonian is expressed as 21] 



^ ^ '^mcr^mm'Cm'cr', (3) 



mm' (7(7' 



where hard wall boundary conditions account for con- 
finement on the L X L lattice with lattice spacing a 
(typically a ~ 3 nm 0). Here cl^^ (Cmo-) is the cre- 
ation (annihilation) operator of an electron at the site 
m = {mx,my). The generalized nearest neighbor hop- 
ping t^^, = (tmm')c<T' accounts for the Rashba couphng 



-tols - itsoc^y (m = m' + e^,) 
-tols + itsoo-x (m m' + Gj^) 



(4) 



through the SO hopping parameter tso = a/2a {Ig 
is the unit 2x2 matrix in the spin space). The di- 
rect correspondence between the continuous Eq. (jSJ and 
the lattice Hamiltonian Eq. Q is established by using 
for the orbital hopping and by select- 
ing the Fermi energy Ep — — 3.8io close to the bottom 
of the band at —AMo to ensure the parabolic energy- 
momentum dispersion. 

At time t' = — oo the 2DEG and the leads are 
not connected, while the left and the right lead are in 



their own thermal equilibrium with the chemical po- 
tentials /ii, and iJLn, respectively, where piL = /^^ + 
eV . The adiabatic switching of the hopping parame- 
ter connecting the leads and the 2DEG generates time 
evolution of the density matrix of the structure [l^ . 
The spin accumulation is obtained as the nonequilib- 
rium statistical avera ge ( ...) (with respect to the den- 
sity matrix at = [ij) of the spin-density operator 

Sm) = f Y.aa' ^ 'y^' (cLijCmCT'), which is expressed via 



the lesser Green function (cJ^jg-CmV') 



i m'm.cr'cr 



0) = ^^IZ.dEG^,^^^,^iE) thereby yielding Eq. ©. 
Here we utilize the fact that the two-time correlation 



(t) ) depends 
it 



function G^^,^^^,{t,t') = 

only on T = t — t' in stationary situations, so 
can be Fourier transformed to energy G^^, (r) = 
A f°° dEG< .{E)e'^-'/^. The matrix G<(£;) 

ZTTh J— oo mm ,crcr \ / v / 

is obtained from the Keldysh equation G^(£^) = 
G''{E)T,<{E)G''{E). Within the effective single-particle 
picture, this equation can be solved exactly in the nonin- 
teracting electron approximation by evaluating the re- 
tarded G-'iE) = Ji; - H - - S£ - Sy-i and 
the advanced G°-{E) = [G^{E)\^ Green function ma- 
trices. In the absence of inelastic scattering, the re- 
tarded self-energies T,1{E - ey/2) and SU^ + eV/2) 
introduced by the left and the right lead jlj], respec- 
tively, determine S<(i;) = -2i[l\rLTs l{E - eV / 2) f l{,E - 
eV/2) + luiHaiE + eV/2)fR{E + eV/2)] [Jlm are the 
Fermi distribution functions of the electrons injected 
from the macroscopic reservoirs through the leads and 
ItsiT,l,r = (S£ — Sf^jj)/2i]. The gauge invariance of 
measurable quantities with respect to the shift of electric 
potential by a constant is satisfied on the proviso that 
S£, depend explicitly on the applied bias voltage V 
while G^{E) has to include the electric potential land- 
scape within the sample (2^. 

Non-equilibrium spin accumulation in clean 
mesoscopic two-probe 2DEG structures. — The one- 
dimensional spatial transverse profile of (S(a; = xq, y)) 
across the macroscopic (200a x 200a) 2DEG sample sup- 
porting steady-state charge current at zero-temperature 
is plotted in Fig. El for the nonlinear eV ^ Ep as well 
as the linear eV <C Ep transport regime. In the linear 
phase-coherent transport regime at low temperatures 
only states at the Fermi energy contribute to (S(r)) in 
Eq. which therefore encodes the information about 
their wave function and oscillates within the sample. 
Figure El also shows that the width of the edge peaks 
of (^'^(r)} is determined by the spin precession length 
Lso — 7rato/2tso [13 (on which spin precesses by an 
angle tt; Lso also plays the role of the DP spin relax- 
ation length 0, in weakly disordered systems) in our 
two-probe device which is in the regime L > Lgo 3> A^?. 

While experiments 0] have searched for the out- 
of-plane spin accumulation exhibiting opposite signs on 
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the lateral edges as the evidence for the spin Hall effect, 
we predict that both {Sz{r)) and {Sx{r)) satisfy this cri- 
terion. The emergence of {Sx{r)) 7^ is the hallmark 
of the mesoscopic spin Hall effect in Rashba wires 
where Hamiltonian Eq. Q gives rise to the SO force op- 
erator [l3 F = 2a'^m*{p x z) ® a^/hP ~ dVconi{y)y/dy 
which deflects spin-t and spin- J, electronic wave pack- 
ets in opposite transverse directions, while spin is at the 
same time precessing since injected |t), | i) (the eigen- 
states of &z ) are not the eigenstates of the Zeeman term 
<T • Br (p) whose BR(p)-field is almost parallel to the 
y-axis 0, 0| . This semiclassical picture 0| heuristi- 
cally explains the symmetry properties of the quantum 
transport induced accumulations in Fig. |21 with respect 
to the bias voltage reversal, {Sz{r))_y = — {Sz{r))y 
vs. {Sx{r))_Y — {Sx{r))Y, which can be tested exper- 
imentally JO]. Also, the a^-dependent transverse SO 
"force" accounts for the difference (S'z(r)}_^ — {Sz{r))^ 
vs. (5,(r))_„ ^-{Sxir))^. Finally, Fig. Ha) confirms 
that (5*2 (r)) spin densities will move away from the lat- 
eral edges upon entering the right lead where a = and 
the transverse SO force deflecting the spins is absent. 

In contrast to the out-of-plane (5*2 (r)) and the in-plane 
longitudinal {Sx{r)) spin Hall accumulation, the trans- 
verse in-plane spin accumulation {Sy{r)) has the same 
sign on both lateral edges. Thus, it cannot originate 
from the spin Hall effect. Instead, it is an analog of the 
magneto-electric effect — where electric field induces spin 
polarization of conduction electrons and nonequilibrium 
magnetization M = /3E — which was argued to occur in 
disordered paramagnetic systems (metals with SO scat- 
tering off impurities [2^ or diffusive Rashba spin-split 
2DEG|l9|,|25j), and has finally been observed experimen- 
tally |26j. It arises due to the combined effect of SO cou- 
pling, inversion asymmetry, and time-reversal symmetry 
breaking by electron scattering (which is necessary due 
to M being t-odd while E is t-even). However, here we 
predict that passage of electric current through ballistic 
sample (E = inside 2DEG) will also induce the trans- 
verse in-plane magnetization, where finite quantum point 
contact conductance of the nanostructure in Fi g. [J (i.e ., 
the corresponding dissipation in the reservoirs [lj,ll5j) 
signifies the time-reversal symmetry breaking. 

The key issue for spintronics applications is to detect 
(at least indirectly the pure spin Hall current fiow- 

ing out of multiterminal structures. When transverse 
leads (e.g., labeled by 2 and 3) are attached to the lat- 
eral edges of 2DEG, the nonequilibrium spin accumula- 
tion will push the pure spin current ^ (J-J —I2) into them, 
as shown in Fig. 0(b). The properties of its spin Hall 
mY = -[7|]^ [J|]^ = -[IIY and spin Warization" 
[/|]^ = [/|]^ components predicted in Ref. 'f?| are exactly 
the same as can be inferred from the profiles of (S(r)} in 
Fig.H 

The effect of disorder and in-plane magnetic field on 
the mesoscopic spin Hall accumulation. — We show in 
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FIG. 3: (color online). The transverse profile {Sz{x — 13a, y)) 
of spin Hall accumulation as a function of spin-independent 
disorder W setting the mean free path £ « 21.5a/VF^ [the 
ballistic limit W = is the transverse profile of (S'z(r)) in 
Fig. H^a)]. The inset shows the effect of the external in- 
plane magnetic field g'/^sBy"* on {Szi^ ~ 13a, y = )) at 
its maximum j/™^" = 26a or minimum t/min = 5a peak values. 



Fig. 121 that disorder, introduced as the random poten- 
tial em e [-W/2,W/2] in Eq. ©, is gradually dimin- 
ishing the amplitude of the original clean limit edge 
peaks of (Szijc)) while leaving their position intact. Thus, 
the magnitude of the peaks alone does not provide sin- 
gle criterion to differentiate between possible underly- 
ing mechanisms of the spin Hall effects 0, ^|. We of- 
fer here yet another experimental tool (in addition to 
{Sx{r))_y = {Sx{r))y test) to identify the mesoscopic 
spin Hall mechanism by demonstrating in the inset of 
Fig.|21that an external in-plane transverse magnetic field 
By^*^, which in strictly 2D systems does not affect or- 
bital degrees of freedom, will not destroy the edge peaks 
of (S'z(r)), except when the Rashba field Br(p) term is 
cancelled by the Zeeman spin splitting term g* ^bB^^ . 
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